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Abstract. Conditions for the establishment of small den- 
sity perturbations in a self-gravitating two component 
fluid mixture are studied using a dynamical system ap- 
proach. 

It is shown that besides the existence of exponentially 
growing and decaying modes, which are present for values 
of the perturbation wave-number k smaller than a critical 
value kj^j , two other, pure oscillatory, modes exist at all 
scales. For k < k^^^ , the growing mode always affects both 
components of the fluid and not only one of them. 

Due to the existence of a resonance between the bary- 
onic and the dark perturbations, it is shown that the onset 
of structure formation in the post recombination epoch is 
substantially enhanced in a narrow scale band around an- 
other critical value fee- For dark matter particles having 
a mass ~ 30 eV, the corresponding critical mass scale for 
the establishment of density perturbations at the time of 
recombination is of the same order of magnitude as the 
galactic one. 
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1. Introduction 

The study of the mechanisms responsible for the forma- 
tion of structures in the universe was started by the pio- 
neer work of Jeans (1902, 1928). Studying the conditions 
under which a cloud of gas (in a static background) be- 
comes gravitationally unstable under its own gravity Jeans 
concluded that perturbations with masses smaller than a 
"critical" value Mj (the Jeans mass) do not grow and 
behave like acoustic waves whereas perturbations with a 
mass M bigger than M ^ grow under the effect of their 
self-gravity, leading to gravitationally bound structures. 
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Since then, the mechanism proposed by Jeans has been 
extensively applied as a criterion of stability in several 
models of galaxy formation. However, by the recombina- 
tion epoch (which is believed to be the one after which 
baryonic perturbations could begin to grow), the criti- 
cal scale predicted by the classical Jeans theory is not at 
all related to the galactic scale. Instead, it is well known 
that just before recombination the critical Jeans mass is 
~ 1O^^~^^M0, (i.e. the mass of rich clusters or even su- 
perclusters of galaxies) , and immediately after recombina- 
tion the Jeans mass drops abruptly, more than 10 orders 
of magnitude, to a value ~ 1O^~^M0, which is related to 
the mass of globular clusters (Weinberg 1972; Zel'dovich 
& Novikov 1983). 

The origin of this discrepancy between the Jeans mass, 
before and after recombination, and the mass of a typical 
galaxy, a few IO^^Mq, has not been, up to now, well un- 
derstood. 

One should emphasize, however, that these values for 
the Jeans mass just before and just after recombination 
epoch are obtained assuming a baryon-dominated uni- 
verse. 

On the other hand, the observation of flat galactic ro- 
tation curves at large galactic radii as well as high galaxy 
velocity dispersions in clusters has led to the missing mass 
problem and to the dark matter conjecture to solve it. 
This fact has, in particular, led to the hypothesis that 
dark matter halos exist around galaxies. 

The existence of dark matter therefore also implies 
that when studying the dynamics and the formation of 
structures in the universe one cannot use the simple Jeans 
stability criterion for a one component gas, but one has 
to study how does the gravitational instability arise in a 
mixture of at least two components. This formulation of 
the problem leads to a different Jeans mass, as well as 
to a better understanding of the origin of galactic mass 
spectrum. This will be the main goal of this paper. 

We must note that, after we finished this work, it was 
pointed out to us by Dr. Varun Sahni that this prob- 
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lem had already been studied, using a different formal- 
ism, by some russian authors like Grishchuk & Zel'dovich 
(1981) or Polyachenko & Fridman (1981). Although hav- 
ing reached some of our conclusions, such studies lack most 
of the features present in this work. In particular they do 
not refer to the existence of the resonance we shall point 
out in this paper. 

Since it is generally accepted that dark matter is con- 
stituted by WIMPs, (Weakly Interacting Massive Parti- 
cles) in this study we shall consider the case for which the 
two components of the cosmic fluid interact only gravita- 
tionally. 

The two fluid components that we shall consider will 
be a baryonic one and a hot dark matter (HDM) one, here- 
after component B and D respectively. The latter being 
assumed to be made of neutrino-like particles. 

Although the numerical results obtained in this paper 
refer to the particular case of HDM, one should point out 
that the formalism here developed is quite general and 
applies to any type of a two-component mixture. 

In Sect. ^ we shall establish the linear coupled equa- 
tions which describe the dynamical system for the den- 
sity perturbations in the two-component fluid we wish to 
study. We shall then develop, in Sect. |[ a qualitative anal- 
ysis of the dynamical system. 

In Sect. ^ we shall study the behaviour, in phase-space, 
of the trajectories representing the general solution of the 
dynamical system. This analytical study will be comple- 
mented with some numerical plots in order to clarify some 
of the qualitative results obtained. 

For numerical purposes we shall consider that the on- 
set of gravitational instability occurs in a flat {il = 1) 
universe at recombination epoch, after the decoupling be- 
tween radiation and baryonic matter, when the radiation 
temperature is ~ 3000 K (Kolb & Turner 1990). We shall 
assume that the density parameter fis has a value in 
agreement with the limits imposed by standard nucleosyn- 
thesis (Walker et al. 1991). Assuming 0.5 < h < 1, where 
h is the Hubble constant in units of 100 Km/s/Mpc, i.e. 
h = Ho/(100 Km/s/Mpc), such a value is - 0.05. 

The neutrino-like particles which constitute the dark 
matter component are assumed to have a mass compati- 
ble with the assumption of being already non-relativistic 
by the recombination epoch. It will be shown later in this 
paper that this value for the mass of the dark matter par- 
ticles seems to provide a good fit for the galactic mass 
spectrum. 

In Sect. I we shall describe the occurence of a reso- 
nance which, in our opinion, is responsible for the forma- 
tion of structure at the typical galactic scales. 

The relation between the existence of the resonance, 
the galactic scale and the mass of dark matter particles is 
discussed in Sect. ||. 

Finally in Sects. ^ and ^ we shall point out the main re- 
sults presented in this paper and outline some open ques- 



tions as well as the research paths we intend to follow in 
the near future. 

2. Dynamical equations for a two component fluid 

The dynamics of a non-relativistic self-gravitating fluid 
can be described by the usual set of hydrodynamical equa- 
tions: continuity equation, Euler equation and Poisson's 
equation. 

In this paper we shall neglect the expansion of the uni- 
verse, since we shall be interested only in the qualitative 
behaviour, at the onset of the process of growth of density 
perturbations in the post recombination epoch, and not 
in the exact dynamics of the process. 

In this case, for a fluid made of two components, D 
and B, the above set of equations becomes: 



dt 



dt 



V-(p«t;«) = 



-^ + {v^.V)v^^ Vp^ -V$ 



dt 



V2$ = 47rG(p„+pJ 



(1) 
(2) 
(3) 

(4) 
(5) 



where $ is the total self-gravitational potential, pj, Vj and 
Pj are respectively the density, velocity and pressure fields 
for fluid elements of component j, (the index j taking the 
values D or B) . 

Using a perturbative analysis, to first order in the per- 
turbation, one can write the dynamical variables as fol- 
lows: 



V — V 5v . 

J JO J 



(6) 



where the index "0" indicates the unperturbed value, and 
6 the first order perturbation. 

It's common to consider the zero order solution to be 
an infinite fluid at rest, {v^ = 0), for which the density and 
pressure are the same everywhere (unperturbed fluid) . We 
must point out, as it is well known, that this static solu- 
tion for the unperturbed state of the above equations is 
not mathematically correct, except for = 0. However, 
since this problem is removed when considering the more 
realistic case of a fluid embedded in the expanding uni- 
verse, (which we shall discuss in a forthcoming paper), 
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one can start perturbing this static configuration in or- 
der to obtain some quafitative information on the fluid 
behaviour. 

We shafl assume that the solutions of the above cou- 
pled system of Eqs. (0) - (||) are square integrable func- 
tions of the spatial variables and therefore can be Fourier 
analysed with respect to these variables. 

To first order in the perturbation, after the Fourier 
analysis in the spatial variables and some manipulation, 
the set of Eqs. (|l|) - (||) leads to the following equations: 



components are only functions of the parameter k and has 
the form: 



A = A{k) 



r 




1 










Vo 














W„- cl P 



1 







(11) 



(7) 



where the double dot stands for second time derivative, k 
is the wave number of the perturbation. For each of the 
components j, the parameter Wj is given by Wj = AnGp^^ 
and Aj is the density contrast Sp. /p. and Cj the adiabatic 

1/2 

sound speed given by — {dp J dp 

The value of the parameter Wj , at the cosmic time of 
recombination, trecj can be related to the density param- 
eter, of component j, by 



(8) 



We shall be particularly interested in the solutions 
which represent the density perturbations in components 
D and B. These solutions are given by the second and 
fourth components, X2(t) and X4{t), of the general solu- 
tion x{t), of the dynamical system (^0|). 

Using methods of qualitative theory of differential 
equations one can analyse the dynamical system and ob- 
tain some qualitative information on the behaviour of the 
perturbations. This will help understanding the stability, 
against density perturbations, of the two-component sys- 
tem. Even in a case like this one where an explicit solution 
exists, this is worth doing. 

Alternatively the autonomous dynamical system can 
be numerically studied provided we specify the initial con- 
dition Xq = a; (to), where to is the cosmic time at which 
the equilibrium state of the fluid is perturbed. 

We shall analyse the system qualitatively and, when- 
ever it proves convenient, to clarify any result, we shall 
illustrate this study with some numerical plots. 



The adiabatic sound speed in component j at recom- 



bination, cj„^, is the sound speed of a monoatomic ideal 3. Qualitative study of the dynamical system 



gas, which is given by: 



3 



,T,. 



1/2 



(9) 



where K^^^ is the Boltzman constant and Tj^^^ and rrij are 
respectively the temperature at recombination and mass 
of the particles of component j. 

One should emphasize, that from the decoupling be- 
tween the neutrino-like particles and radiation until the 
epoch at which these particles become non-relativistic, 
(which occurs for a red-shift ^ 6 x 10'*, Doroshke- 
vich et al. 1988), their temperature T^i is given by T^ = 
(4/11)^/^ T^, where T^ is the radiation temperature. Since 
then, it decays with a^^, where a is the scale factor of 
the universe (Gao & Ruffini 1981; Padmanabhan 1993). 
One can therefore estimate that, at the recombination 
epoch, such a dark component D will have a temperature 
Tz3„, ~ 40 K 

The above system (0) of two linear differential equa- 
tions of second order can be transformed into the following 
autonomous dynamical system of four first order differen- 
tial equations which one can write in a matrix form as: 



The above matrix A{k) has 4 different eigenvalues for all 
values of k except for fc = and for a critical value k = k^,, , 
which shall be defined below. Therefore for fc > and 
k 7^ k^,j there are 4 linearly independent eigenvectors. In 
this case the general solution x{t) of the system (10) is 
weU known (Arnold 1973). 

If all the eigenvalues are real ones, that solution can 
be written in the form of a linear expansion of particular 
solutions c^J*^ -, i.e.: 



x[t) 



4 



:e^^*^, 



(12) 



X = Ax 



(10) 



where x = (xi, 12, is, 0:4)'^ = (A^ , , , A^)'^, the 
superscript T standing for transpose. A is a matrix whose 



where aj are real-valued functions of the the parameter 
fc, which are determined by the initial conditions, and 
and ^j are respectively the eigenvalues of the matrix A{k) 
and its associated eigenvectors. 

If some of the eigenvalues are complex conjugates, one 
can easily see that the real and imaginary parts of the 
complex function c^J * ^j are two real- valued particular so- 
lutions of the dynamical system associated with the com- 
plex conjugate pair of eigenvalues Xj and Xj. 

In order to obtain some qualitative insight about the 
behaviour of the perturbations, one must find the criti- 
cal or equilibrium points of the system (|lO|) and calculate 
the eigenvalues of matrix A{k) as well as their associated 
eigenvectors. 
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3.1. Critical points, eigenvalues and eigenvectors 

The critical points of the dynamical system (|l^ are ob- 
tained solving the equation 



a; = 0. 



(13) 



The only critical point which does not depend on the 
parameter k is the origin, x = (0, 0, 0, 0)"'^. 

The parameter-dependent critical points are given by 
the solutions of the equation det(A) = 0. In the case when 
the parameter k has a value given by: 



(14) 



any point x = (0, Xj, 0, x^)""" of the (x2,X4) plane, is a 
critical point of the above system. 

In Eq. , fc^ and kj^ are respectively the Jeans wave 
numbers of components B and D when taken separately, 
and we shall call k^^j the Jeans wave number for the mix- 
ture of the two components. 

The eigenvalues of the matrix A{k) are the roots of the 
characteristic polynomial, obtained solving the equation 



det(A-AJ) = 



(15) 



where / is the identity matrix. 

The 4 eigenvalues of A(fc), which we shall represent by 
A 1 , A 2 , A3 and A 4 , are therefore given by: 



Al 



A2 = ^^^z + VT^+iPff 



(16) 



,A3 = -A4 = -^^Jf-^P + ^g 
where / and g are the following functions of k: 
f(k)^W,+W,-k^cl+cl) (17) 



gik)^W,cl+W,cl-k' 



(18) 



The eigenvectors of the matrix A{k), associated to 
the eigenvalues A^, must satisfy the condition 



and can be written in the following form: 

= (/3iAi, A, Al, 1)T 
I, - (/32A2, /52, A2, 1)T 
= (/33A3, [33, A3, 1)^ 

,^4 = iPi A4, P4, A4, 1) 

where Pj are also functions of k and are given by: 

>i = P2 ^ 2w^[h + ^WTmTnr 

and h = h[k) is given by: 

h{k) = w^-w^+e (4-4) 



(19) 



(20) 



(21) 



(22) 



Table 1. The signs of real and imaginary parts, (5R and 9 re- 
spectively), of the four eigenvalues of matrix A{k) as a function 
of parameter k. Note that S$(A3) = — £s(A4) 7^ for all k ^ Q. 
(See text for discussion.) 
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3.2. Qualitative results 

Qualitative information about the behaviour of the dy- 
namical system can be obtained studying the signs of the 
real and imaginary parts of the eigenvalues Xj . 

The only zero of the eigenvalues A 1 and A 2 is obtained 
when k is the Jeans wave-number of the mixture, (fc = 

Looking at Eqs. (|l|) - (|l|) one can see that for fc < fc,, 
the function ^(fc) is positive and therefore Ai and A 2 are 
real nonzero eigenvalues, (Ai > 0, A2 < 0). For fc > fc^,,,, 
Al and A 2 are two pure imaginary conjugate eigenvalues. 

The eigenvalues A3 and A 4 do not have zeros, except 
for the asymptotic value fc = 0. Since / — \J P + Ak'^g is 
always negative, for all other values of fc A3 and A 4 are 
pure imaginary conjugate eigenvalues. 

Physically, this means that one has two paired wave 
modes (which correspond to acoustic oscillations), for all 
values of the wave number fc and two other modes which 
are acoustic for fc > fc^^- and become one growing and one 
decaying mode for fc < fc^^ . 

These results are summarized in Table |^. 

One can therefore distinguish two, very distinct, sit- 
uations concerning the four-dimensional phase-space ac- 
cording to the values of the parameter fc, i.e.: 

1. For fc > fcjj- the four eigenvalues Ai, A2, A3 and A4, 
are pure imaginary and therefore the phase-space is 
a center space, which is spawned by the eigenvectors 
associated to those eigenvalues. 

Solutions in a center subspace oscillate at constant am- 
plitude, and therefore density perturbations with wave 
number fc > k^^^ can't grow. This case corresponds 
to the existence of two pairs of acoustic wave modes 
propagating in the fluid with different phase veloci- 
ties. In the limit of small wavelengths, it is easy to see 
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that these modes have the same dispersion relations 
and compositions as two completely uncoupled acous- 
tic waves propagating in pure D and B fluids. This 
is due to the fact that the role played by their self 
gravity, which is the only mechanism coupling the two 
components, is not important at small scales. 
2. For k < k^, (see Table ||), the eigenvalue Ai is real 
and positive, A2 is real and negative and A3 and A4 
are pure imaginary (conjugates) eigenvalues. In this 
case the phase-space is, therefore, the direct product 
of three distinct subspaces: 

(a) An unstable subspace of dimension one spawned by 
the eigenvector associated to the eigenvalue with 
positive real part, (Ai). 

Solutions in an unstable subspace always grow. 
Since Ai is a real time-independent eigenvalue the 
growth of these particular solutions is exponential 
without oscillations. 

(b) A stable subspace also of dimension one spawned by 
the eigenvector associated to the eigenvalue with 
negative real part, (A2). 

Solutions in a stable subspace always decay. In our 
case the eigenvalue associated to this subspace, A2, 
is also a real valued one and therefore the decay 
of these particular solutions is exponential without 
oscillations. 

(c) A center subspace which for this range of k is of 
dimension two and is spawned by the eigenvectors 
associated to the complex eigenvalues with null real 
part, (A3 and A4). Solutions in a center subspace 
oscilate. 

This case corresponds to the existence of a pair of 
acoustic modes plus one growing and one decaying 
mode. 

The analysis of the composition of these modes shows 
that all of them envolve both components D and B, 
and are completely coupled through gravity, which, at 
these large scales plays an important role. 

One can then conclude that growth of density per- 
turbations is possible only for /c < fc„. Note that for all 
k < kj^ this growth affects the two components (i.e. the 
mixture) and not only one of them. This is the reason why 
we called k^^ the Jeans wave number of the mixture. The 
corresponding physical critical scale length below which 
no structure could develop is ^„ = 2-1: /k^^. 

It is interesting to note that although this scale length 
is related to the Jeans lengths, 1^ and l^ , of the two com- 
ponents when taken separately it does not coincide, in the 
general case, with none of them. In fact from relation ( p^ 

one can see that 1^^ = (^^■^ + Z^-^) 

Since we are mainly interested in structure formation 
we shall consider, in what follows, only the case when k < 
fc„, (i.e. I > l^,). 



4. Behaviour of the general solution 

In order to study the behaviour with scale of the dynami- 
cal system at the onset of the structure formation process, 
one needs to analyse the /c-dependence of the eigenvalues 
Aj, the eigenvectors and the coefficients aj, which ap- 
pear in Eq. (|l2|). 

Since A2 — — Ai and A4 = —A3 we shall only study the 
dependence with scale of the real eigenvalue Ai and of the 
imaginary part of A3 . 

The eigenvalue Ai, being real and positive, is the one 
related with the growth rate of the perturbations, is a 
growing function of scale. However Ai does not grow at 
the same rate for all scales. In fact there are two particu- 
lar scales at which a significant change in its growth rate 
occurs. The first one is the Jeans scale-length of the mix- 
ture, Ij^j, where, after changing from an imaginary to a 
real value, Ai increases very fast and then slows down, re- 
maining approximately constant until the second peculiar 
scale is reached, which occurs near the Jeans scale-length 
of the dark matter component, where Ai grows faster 
again. Beyond this region the growth of Ai slows down 
again and approaches the asymptotic value + as 
/ ^ 00; see Eqs. (^ and (|l^). 

The eigenvalue A3 is a decaying function of scale and 
the only change in its behaviour occurs at a scale which 
coincides with the second peculiar scale for Ai, (i.e. for a 
scale near l^), where a faster decrease than for any other 
region of the scale range, is observed. For greater scales 
it approaches zero as / 00. Since A3 is the frequency 
associated to the acoustic modes one can see, from its de- 
pendence with scale, that although oscillations are present 
in all scales, for large scales acoustic waves have their fre- 
quency approaching zero, i.e. their period tends to infinity. 
See Fig. 0. 

Additional information about the system behaviour 
can be obtained studying the trajectories in phase-space 
corresponding to the general solution x{t); unfortunately, 
since the phase-space is of dimension four, they are hard to 
visualise. However, since, as referred above, we are partic- 
ularly interested in the qualitative behaviour of the second 
and fourth components of the general solution x{t), i.e. the 
behaviour of the density perturbations A^(t) and A^(t), 
we shall infer some properties of the dynamical evolution 
of the density perturbations analysing the behaviour of 
the projection, on the plane (a;4,X2) = (A^, A^-,), of the 
phase-space four-dimensional trajectories. 

Looking at the components of the eigenvectors, given 
by the relations ( |l6|) - (22), one can see that the projection, 
on the plane (0:4, 2:2), of the stable subspace coincides with 
that of the unstable subspace (since j3i = (32)- 

The common projection of these subspaces is the 
straight line given by the equation: 



A„ =/3iA, 



(23) 
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The center subspace has a projection on the same plane 
(Ag, A^) which is also a straight line, defined by 



(24) 



These equations, ( p3|) and (pij), for the projection of 
the stable/unstable and center subspaces can give us the 
relative composition (relative percentage of components 
B and D), in the growing and decaying modes, Eq. (p3|), 
as well as in the acoustic modes, Eq. (|24|). 

From relations (|2^), one can see that since + 
AW^W^ is always positive, (3j are, therefore, real- valued 
functions for all k > 0, with /3i > and /Js < 0. The 
condition /3i > means that both components partici- 
pate positively in the collapse of the mixture as well as in 
the decaying mode, which being transient will not interest 
us. On the other hand < means that, in the acoustic 
modes, the two components oscillate with opposite phases; 
when the density in one of them is growing, the density in 
the other is decreasing. 

The direct product of the stable and unstable sub- 
spaces forms a two dimensional subspace which is usually 
called a saddle subspace and solutions in such a subspace 
are simultaneously pulled towards the equilibrium point 
along the stable axis and pushed away from it along the 
unstable axis. The result is that when the time t varies 
from — oo to +oo the orbits representing these particular 
solutions are initially attracted by the critical point reach- 
ing a closest approach to the origin at a time td after which 
they are repelled away. 

Choosing initial conditions such that xi = X3 = 0, i.e. 
Xq lies on the plane (2:4, 0:2), and given the symmetry, in 
first and third components, between the eigenvectors ^ i 
and oiic concludes that the closest approach from the 
origin of these orbits, is attained at initial time to, (i.e. 
tci = to)- Therefore, in this case, for t > to the system is 
monotonically driven away from the origin with time. It is 
worth noting that this condition (xq lying on the [x^, X2) 
plane), corresponds to the reasonable assumption that the 
mixture starts to collapse from rest. 

Moreover, since the effect of the center subspace is to 
force the orbit to oscillate, one shall have a projected path 
which oscillates around the projection of the saddle sub- 
space, i.e. around the straight line given by Eq. (p3[), see 
Fig. |. 

According to the initial conditions, A „ and A„ , one 
can see that, if A^^ > /Js A^^, the system will collapse 
to form a massive structure. On the other hand, if A„ < 
/?3 A it will evolve to form a void, as represented in Fig. 

Since (3i and (3^ depend on the perturbation's wave- 
number k, it is clear that the direction of the projection 
of the above mentioned subspaces also changes with scale, 
and therefore the relative growth of the perturbations in 
the two components, as well as the composition of the 
collapsing and oscillating modes, is also scale dependent. 



The change in the relative growth of perturbations 
with scale can be seen analysing the behaviour oi Pi with 
wave number k, or equivalently with scale length I. 

The straight lines defined by Eqs. ( p3| ) and ( [24[ ) make 
angles with the direction of the a;4-axis, which are given 
respectively by 9i — arctan(/3i) and 9^ = arctan(/33). 

If one assumes that < (case of hot dark matter) , 
the derivative 



dl W„ 



1 



(25) 



is always positive and since (3i{l^j) = '^b/'^d ^'^"^ 
limj^oo /?! = 1 one concludes that di is a monotonic func- 
tion of I, growing from arctan(c^ /c^ ) imtil the asymptotic 
value limj^oo 61 — 45°. 

By a similar study oi (3^ one concludes that 



dp 
"dl 



' ^ 4^2 



1 - 




(26) 



is also always positive and therefore 9^ grows with scale 

w c 



from 93(1 j^j) — arctan 



until the asymptotic 



arctan 



(-Sf) 



This is illustrated in 



value lim /^oo 6 
Fig. I 

It is important to note, however, that the rates of 
change with scale of both jSi and /^a are not constant. 
The same is obviously true for 9i and ^3 as shown in 
the plot of Fig. ^ In fact from this plot, or alternatively 
analysing expressions ( p5| ) and (^6|), one can notice that 
for Ij^j < I <^ (i.e. fcj,j > k ^ k^), Pi and P^ are slowly 
increasing functions, while when I reaches the order of 
they become much steeper functions of I. For values of I 
above this scale the rate of change oi Pi and P^ slows 
down again as Pi and P3 approach the asymptotic values 
1 and — Wg /W^ respectively. 

This means that the composition of the mixture in the 
collapsing mode changes substantially when the pertur- 
bation scale is of the order of , starting to involve, at 
this scale, a much greater part of the dark material which, 
for smaller scales, was mainly oscillating in the acoustic 
mode. 

Further on, it will be shown that there is a critical 
scale L where a resonance between the two components 
occurfjd. It is for this particular scale, and larger ones, 
that the acoustic modes tend to be devoided in favour of 
the collapse mode. 



^ This scale Ic is in the particular case we are considering of 
the order of as shall be seen later. Moreover it will be shown 
that, if the mass of hot dark matter particles is m^-, ~ 30 eV, 
the scale Ic has a corresponding mass of the same order as that 
of typical galaxies. 
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As a matter of fact, in the limit of very large scales 
(fc 0) one has: 



lim ^ = 1 



(27) 



where A^^^ and A^^ are the density contrast of the dark 
and baryonic components participating in the collapse 
mode. 

Denoting by and (5^^ the density perturbations of 
the dark and baryonic components in the collapse mode, 
one has: 



lim 



(28) 



and for the acoustic mode, by a similar reasoning, one 
obtains 



lim P^ = -l 



(29) 



where S^^ and (5^^ are the density perturbations in the 
dark and baryonic components in the acoustic modes. 

These results mean that, for very large scales, the com- 
position of the perturbed region, in the collapsing mode, 
is identical to the one in the unperturbed state, as shown 



by Eq. (28). On the other hand the perturbations in both 



components, oscillating in the acoustic modes, cancel each 
other, since they oscillate with opposite phase and the 
same amplitude. 

5. The structure formation resonance 

The coefficients aj in Eq. ( p^ ) can be obtained by spec- 
ifying the initial conditions Xi{to) — Xi^. Substituting in 
the general solution given by Eq. (p^), one obtains: 



2^ 



(30) 



where iifj), is the i*'* component of eigenvector j . 

Solving the above linear system of equations along with 
the already mentioned assumption Xi^^ = ccsg = one 
obtains the coefficients aj , which are given by: 



"1-2 /3i-/33 ^ 



_ fHo 1-Qo03 Alto 
"2-2 /3i-/33 ^ 



as 



= ^2o COs(zA3to) 



(31) 



"4 = X2a ^j^^ sin {iX^to) 



where Qq is the ratio between the density contrast in 
the two components at i — to, (i.e. Qq = X4„/x2„ = 
A„ /A„ ), and i = 



Substituting these expressions in the second and fourth 
components of the general solution ([l^), one obtains after 
some straightforward manipulation: 



A^(t) = X2{t) 

= X2o [Ci (e^^" + e-^'-) + (2 cos(zA3t)] (32) 

and 

A^(r) = X4{t) 

= X4„ [Ca (e^^" + e-^^") + C4 cos(iA3T)] (33) 

where we have used a new time variable, r = t — to, by 
translating the time origin to the beginning of the process. 
The new coefficients Q, which are the mode amplitudes 
for the perturbations in components D and B, are given 
by the following expressions 



Ci 



2 /3i-/33 



^ _ 1 Qn^-f>3 
^3-2 01-03 



(34) 



C4 



/3i-C 



01-03 



For simplicity it shall be assumed, throughout this sec- 
tion, that Qo is positive, although massive structures can 
also originate even in the cases when it is negative, (see 
the previous section, Fig. ^). 

One must also point out that if one allows xi^ ^ 
and/or x^g ^ the expressions ( ^l| ) for aj, and conse- 
quently the expressions ( ^ ) for Q, become more com- 
plicated, and therefore dificult to be studied analytically. 
The general analytical study of the mode amplitudes, in 
the case when all the Xi^ ^ 0, is under development and 
we intend to publish it in a forthcoming paper. However 
preliminary numerical results seem to indicate that no sig- 
nificant change occurs in the behaviour of the mode am- 
plitudes. 

Looking at the behaviour of the amplitudes Cj, one 
should first of all point out that, from Eqs. (p^), one can 
see that the pairs (Ci , C2) and (Ca , C4) are related by: 



2Cl+C2 = 2C3 + C4 = l 



(35) 



This means that a maximum of Ci corresponds to a 
minimum of ^2 and a maximum of Ca corresponds to a 
minimum of (^4 or inversely. 

This can be seen in Fig. ^ where the numerial plot 
illustrates the dependence of ^3 on scale, and shows its 
resonant behaviour for I ~ Ic, and in Fig. |6[ where the 
sharp minimum occurring for the coefficient (^4 is observed 
at the same value of the scale length. In the numerical plot 
of Fig. 1^ we also illustrate the behaviour of the coefficients 
Ci and (2- 
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Physically, this means that the occurrence of a reso- 
nance (presence of maxima in or C3), is related with 
the fact that at the resonant scale the ratio between the 
number of particles which populate the collapse and de- 
caying modes and the number of particles of the same 
component which populate the acoustic mode grows sub- 
stantially when compared with the same ratio in nearby 
scales. Therefore one can say that most of the resonant 
material is collapsing and little is left oscillating, with the 
oscillating mode helping to devoid the background to make 
a greater percentage of the particles participate in the col- 
lapse. 

The occurrence of such a resonance can be interpreted 
as a consequence of the following: 

As pointed out by Grishchuk & Zel'dovich (1981), the 
dynamical system (^0[) represents two gravitationally cou- 
pled oscillators. If one takes the two components B and 
D separately, the dispersion relations, for perturbations in 
those separate components, give us the natural frequen- 
cies of oscillation, i.e. the ones the components would have 
if uncoupled, which are given by 



1/2 



and 



OJr 



(cl 



>l/2 



(36) 



(37) 



where and are the natural frequencies for compo- 
nents B and D respectively. 

From Eqs. (^6|) and (|3^) one can see that the two nat- 
ural frequencies coincide for a wave number given by 



1/2 



(38) 



This means that kc is the characteristic wave number 
for which perturbations in the two components, taken sep- 
arately, would have the same collapse time scales. Thus it 
is reasonable to expect a resonance to occur, at least in 
some of the mode amplitudes of the perturbations, for a 
value of k near k^ . In fact it will be shown below that such 
a resonance does indeed occur. 

The scale length corresponding to k^ is given by 



1/2 



il-{wjw^)il 



1/2 



(39) 



In the case when ^ and <C one can see 
that If. is of the order of , as noted in the last section. 

This resonance has physical meaning only if the value 
of kc is real and positive. Looking at Eq. (|3|) one can see 
that this can only happen in the case when the gravitation- 
ally dominant component, i.e. the one with a greater value 
of W , has also a greater sound speed, otherwise the scale 
for which the two components have equal natural frequen- 
cies will be a physical meaningless imaginary one. There- 
fore, in what follows, it will be assumed that W^^ > W„ 



and Cj-, > Cg. This is a necessary condition for the exis- 
tence of a resonance in the system. When applied to our 
universe, for which the observational evidences seem to in- 
dicate to be dominated by dark matter, (i.e. > Wg), 
this means that such a resonance cannot occur in a cold 
dark matter (CDM) dominated universe. 

As can be seen in Eqs. (|3^), the peculiar behaviour of 
(3 1 and /Js for I ~ Ic mentioned in the preceding section, 
is responsible for the fact that, from all the perturbations 
which can start to grow at the recombination epoch, those 
with masses of the order of Mc are the ones which are go- 
ing to present, at the onset of the gravitational instability, 
the biggest amplitude, in the collapse mode. In this sense 
the scale Ic is a privileged one. 

However one must be careful in interpreting this re- 
sult since it is valid only for the onset of the instability 
and cannot be extrapolated in time. To have a correct 
idea about the dynamical evolution of the perturbations 
one has to solve a more complicated problem, i.e. we must 
include, in the equations, the terms due to the expan- 
sion of the universe, which allows us to obtain the correct 
growth rate of the perturbations. On the other hand an- 
other dificulty arises because we solve the equations only 
for one Fourier component of the general solution which 
separates a perturbation into distinct wavelengths that, 
in principle, have distinct growth rates. Therefore at any 
epoch the density contrast in a perturbed region of scale 
^ l/k must be integrated in all the region, i.e. for all k in 
the range < A: < (see Padmanabhan 1993). 

In particular, before a careful analysis of these two 
points is made, one cannot know which perturbations are 
the first to reach the non linear regime. This is a very 
important question that we shall address in a forthcoming 
paper (in preparation). 

Since the growing mode is the one that leads to struc- 
ture formation we shall be particularly interested in the 
behaviour of the coefficients (i and C3 which are the ones 
associated to that mode. 

We shall therefore proceed by analysing the functions 
and ^3. Looking at their derivatives with respect to k, 

one can see that their maxima occur for values of k given 

by: 



for Cl, and 



2W„ 



1/2 



(40) 



2W^ 



-,1/2 



1 (C?, -C2) 



(41) 



for C3, where q = 1/Qa = A^^^A^^^. 
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Substituting Eq. ( p^ ) into Eq. ( pO| ) one can see that the 
existence of a maximum for in the range < k < k-^, 
can only occur if q satisfies the following condition: 



1 A 



(42) 



On the other hand, substituting Eq. into Eq. (|l|) 
one can see that the existence of a maximum for (^3 in the 
same range, can only occur if: 




It is reasonable to assume that the composition of the 
initial perturbation is the same as the composition of the 
unperturbed state (i. e. q = W^/Wg) . In this case, it 
is easy to see that ^1 will never have a maximum in the 
desired range, while ^3 can have one, only if the following 
condition is satisfied: 



(44) 



From this analysis one concludes that, in a hot dark 
matter dominated-universe, the resonance will only affect 
the baryonic component. We shall, therefore, from now on 
study only the behaviour of the mode amplitude (^3 which 
is the one associated to the growing mode in the baryonic 
component. 

Substituting the value of k {(f"^) given by Eq. (|l]) 
into Eq. (^) , one obtains for the function C3 its maximum 
value C^^^, given by: 



^max 
^3 




(45) 



which depends only on the unperturbed densities of the 
components and on the initial conditions. 

In the case of a universe dominated by dark matter 
{W„/Wg ^ 1) and for the initial conditions specified 

\ 3/2 

) 

From this expression it is clear that the height of the 
resonant peak increases as the parameter q, (in this case of 
the order of /W^ ) , increases. This means that the more 
the Universe is dominated by dark matter, the more the 
effect of the structure formation resonance in the baryonic 
material will be important. 

In Fig. 1^ the amplitude ^3 , of the growing mode in the 
baryonic component, is ploted as a function of the mass 
scale of the perturbation, for q = 19, which corresponds 
to assume that q = W^/W^ = Hd/^b, using for fls and 
Qjj the values 0.05 and 0.95 respectively. This is to be 
compared with the plots in Fig. ^ where other values of 
the parameter q were used. 



The importance of the resonance can be also measured 
by the width of its band: a larger resonant band corre- 
sponds to a weaker resonance (Feynman 1966). We shall 
define the resonant band width 7 as the width of the curve 
of C3 at half the maximum of its height. 

The two values of k (fci.2), which define the resonant 
band width 7 = |fci — fc2|, must therefore satisfy the fol- 
lowing condition: 



(43) are given by 



kia 



where 



± 



(%/2-4C; 



max 
3 



J d 



1/2 



(46) 



(47) 



(48) 



and 



q^ + 2 



{V2-1) 



1/2 



(49) 



For given sound speeds of the two components one 
can see that this width is very sensitive to the product 
WgWj-,. On the other hand, the further apart the sound 
speeds and c^^ in the two components are, the smaller 
the resonant band width will be, indicating a stronger and 
sharper resonance. 

Another feature one should point out is the behaviour 
of ^3 to the left and to the right of the resonant scale. 

Since there is no growing mode for k > /c^j , we shall 
be interested in studying the behaviour of (^3 for the limits 
k ^ and k k^., . 

The value of C3 in the limit fc 0, depends on the 
values of Wg and and is given by: 



limC3^^^^- 



(50) 



whereas the value of (^3 for k — k^ depends on both the 
sound speeds and on the densities of the two components, 
and is given by: 



C3(fc..) = 2 



n J s LL 



+ 



w 

B 



(51) 



Assuming, as above, that q ~ Wj^/W^ , the asymptotic 
value limi:^o C3 ^ 1/2, as can be seen from Eq. (^0|), 
while, if one allows q to be much smaller than 1, then 
limfe^o C3 < 1/2. 
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On the other hand, in the case W^/W^ c^/c^, and 
for reasonable values of q, the value of CsC^aj ) always of 
the order of 1/2. 

One can therefore conclude that the curve has a 
different behaviour on the far left and on the far right of 
the resonant scale. For small scales the value of does 
not vary very much with initial conditions, while for great 
scales that value is very sensitive to initial conditions. The 
smaller the initial perturbation in the dark component 
when compared with the initial perturbation in the bary- 
onic component, the smaller will be the value of the mode 
amplitude in the growing mode of the baryonic compo- 
nent. For scales larger than the resonant one, favouring, 
in consequence, the onset of perturbations at the resonant 
scale or smaller (relative to greater ones). The inverse oc- 
curs if the perturbation in the dark component increases 
relative to the perturbation in the baryonic component, 
i.e. in this case the favoured scales for the onset of pertur- 
bations are the resonant one or greater (see FigH). 



6. The resonance, the galactic scale and the mass 
of dark matter particles 

We believe that the resonance described in Sect. ^ is re- 
sponsible for setting up the growth, at the recombination 
epoch, of massive structures at the galactic scale. The res- 
onance therefore provides the mechanism to select out this 
particular scale, among all possible scales which could un- 
dergo the gravitational collapse. In other words this reso- 
nance provides an explanation why do galaxies, with the 
typical masses one observes, form around the recombina- 
tion epoch, even if the Jeans mass for baryons, and even 
the Jeans mass for the mixture, either before or just after 
this epoch, is either much larger or smaller than typical 
galactic masses. 

The typical galactic masses should, according to this 
work, correspond to the mass of the perturbations, which, 
at the recombination epoch, would have a scale near the 
above mentioned resonant scale. One would therefore have 
no need to assume that the mass spectrum of perturba- 
tions at that time was already populated in a way similar 
to the one we observe now. Instead, one can say that the 
typical galactic mass scale would simply result from the 
scale selected by this resonance. 

Of course we need a better and deeper understanding 
of the physics of the resonant mechanism and how it can 
affect the power spectrum of the density perturbations, 
not only at recombination epoch but also in the period 
following it. Therefore one has to go from separate Fourier 
amplitudes to quantities integrated over the wavelength 
scale. Since one is interested in knowing how the spectrum 
of inhomogeneities will transform after recombination, one 
must follow the dynamics of several Fourier components 
of the density perturbations from recombination until the 
epoch when the perturbations become non linear. 



Note that in the matter dominated period, before re- 
combination, no resonance can exist. This is because, due 
to the coupling between baryonic matter and radiation in 
that period, the component with greater sound speed is 
the baryonic one. On the other hand, in a flat universe 
like the one we have assumed, the baryons are not gravi- 
tationally dominant. Therefore one can see from Eq. (38) 
that the resonance cannot occur in that period. 

Thus, in order to have a detailed understanding of 
galaxy formation one only needs to study the dynamics of 
the Fourier components after recombination epoch. Before 
recombination one can assume that the power spectrum 
was any reasonable one predicted by another model, like, 
for instance, some type of inflationary model. 

Since this is a very important point which must be 
carefully analysed. Since it lies out of the main goal of 
this paper, we intend to treat it thoroughly in a separate 
paper. 

One can relate the scale of the resonance with the mass 
of the dark matter particle. The scale at which the reso- 
nance peak occurs depends on , see Eq. ( ^l|) , which in 
turn is a function of the value of the mass of dark matter 
particles . From Eq. ( ^ ) we have already concluded 
that in order for the resonance to occur, the dark matter 
must necessarily be hot. 

The total mass M inside a spherical region of radius 
Z/2 = tt/A; is given by 



fc3' 



(52) 



is the total 



where k is the wave number and p = / 
mean density in that region. 

In order to obtain bounds for the mass of dark matter 
particles we shall assume that galaxies form by means of 
the above described resonance and therefore have mass 
scales which correspond to the resonant scale for C3. The 
range of galactic masses which one observes in the universe 
will therefore imply a specific range for the mass of hot 
dark matter particles. 

Substituting in Eq. ( |52| ) the value of the wave number 
k by its resonant value k{C,™^'^)^ given by Eq. (pl|), one 
can establish the following relation between M (in solar 
masses) and (in eV): 



210 TT^/^ 



1/2^1/3 



(53) 



Pd 



In deriving this approximate equationa, we have as- 
sumed that ^ and that c^^ ^ c^. The sound 
speeds c. at recombination, appearing in Eq. (^ij) are 
given by Eq. (i). 



^ Although an exact, but more complicated, expression can 
be derived for the mass of dark matter particles, the use of 
these assumptions, which are valid for the case of HDM, allows 
us to derive a simple but accurate expression for . 
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One should emphasize that from the decouphng, be- 
tween the neutrino-Uke particles and radiation, until the 
epoch at which these particles become non-relativistic, 
which for ~ 30 eV occurs at a red-shift '--^ 6 x 10'', 
(Doroshkevich et al. 1988), their temperature is given 
by Td ~ (4/11)^^'^ T-y, where Tj is the radiation temper- 
ature. Since then, it decays with a~^, where a is the scale 
factor of the universe (Gao & Ruffini 1981; Padmanab- 
han 1993). One can therefore estimate that, at recombi- 
nation, such a dark component D will have a temperature 
Tr._ ~ 40 K. 

In Fig. ^ is plotted the scale M at which the resonant 
peak occurs (which we have assumed to be the galactic 
total mass) versus the mass of hot dark matter particles 

Assuming that a typical galaxy formed at recombina- 
tion epoch has a baryonic mass of a few lO'^^ M© one can 
estimate from the plot in Fig. || that the mass of the dark 
matter particles is of the order of 30 eV. It is remarkable 
that this value agrees quite well with the cosmological up- 
per bounds for neutrino's masses (Peebles 1993) as well as 
the values predicted by Sciama's Decaying Dark Matter 
Hypothesis (Sciama 1990a; Sciama 1990b). 

7. Conclusions 

It is clear from the results obtained in this work that it 
is incorrect to use the Jeans criterion for a 1-component 
fluid to a 2-component mixture, for two main reasons: 

1. The Jeans scale for the mixture does not correspond 
to any of the Jeans scales for the components taken 
separately. 

2. Moreover, in a 2-component mixture there are always 
present, at all scales, oscillations in the two compo- 
nents, which correspond to acoustic waves. The same 
is not true for a one component system. 

However the main result of this work is, as we have 
shown, that, in a two-component fluid, a resonance be- 
tween the two components of the fluid must always occur 
provided that the component with a greater density is the 
one with a greater sound speed, see Eq. (|38|). One can 
therefore put some constraints on the nature of the dark 
matter if galaxies are to be formed by such a resonant 
effect. 

In fact, assuming that the density parameter of the 
universe is = 1 and adopting a value for the density of 
the baryonic component in agreement with standard light- 
element nucleosynthesis, one is led to the conclusion that 
the universe is dominated by dark matter, i.e. > . 
In such a case, if > (case of CDM), one concludes 
from Eq. (^8|) that no resonance between the two compo- 
nents is possible and therefore, the galactic scale would 
not be a preferred one, as it seems to be. One is there- 
fore led to believe that most, if not all, of the dark matter 
present in the universe is hot, i.e. it is constituted by light 



neutrino-like particles as implied by the condition for the 
existence of a strong resonance, i.e. c„ ^ Cg, (see Eq. (47) 
and comments following it). 

This resonance is linked to the fact that at the critical 
scale Ic, there is an enhancement of the number of par- 
ticles in both components, (but particularly in the bary- 
onic one), which participate in the collapse mode, instead 
of participating in the acoustic modes as it happens for 
other scales. This effect is only significant in a narrow 
resonant band in the length, or mass, scale of the pertur- 
bations. This, we believe, is the reason for the formation 
of galactic structures at the time of recombination with 
their characteristic range of mass scales. 

This effect occurs at the galactic scales provided that 
the dark matter component is made of neutrino-like par- 
ticles with masses of the order of 30 eV. This value for 
the mass of the neutrino-like particles is below the upper 
bounds imposed by cosmological constraints on the neu- 
trino's mass (Peebles 1993), assuming a flat universe, and 
remarkably close to the value predicted by Sciama's De- 
caying Dark Matter Hypothesis (Sciama 1990a; Sciama 
1990b). We must also note that, from the numerical re- 
sults obtained in the preceding sections, this mass scale is 
of the order of lO'^'^M© which gives a baryonic mass, in 
the same scalelength, of the order of -^W^^Mq, which 

lies in the range of the masses of typical galaxies (a few 
lO^Mo). 

One should point out, however, that the resonant scale 
Ic is particularly sensitive to the mass of the particles con- 
stituting the dark component, see Fig. |[ 

8. Future work 

This study has left a lot of open questions, some of which 
are already under investigation. Below we briefly comment 
on some of them. 

First: One should point out that this study assumes 
galaxies to have started their formation during recombi- 
nation. However, this can well be not the case and they 
could have started their formation over a wider time in- 
terval starting from recombination. In that case, the un- 
perturbed densities pj and the sound speeds Cj in the two 
components would vary during such a period, allowing the 
scale of the resonance to be shifted. This might in turn 
provide an explanation for the observed wide mass range 
for galactic structures. This is a problem which we intend 
to address in a future paper (in preparation), where the 
effects due to the expansion of the universe are taken into 
account. 

Second: Although the analysis of possible effects of this 
resonant mechanism on the power spectrum is out of the 
main purpose of this paper, the existence of the above 
mentioned resonance, or at least the change in behaviour 
of the mode amplitudes near the resonant scale is, on its 
own, a very important result. The occurrence of this reso- 
nance can have, in principle, important implications in the 
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process of galaxy formation, and, in particular, transform Zcl'dovich Ya. B., Novikov I. D. 1983, Rolativistic Astrophysics 
the power spectrum of density perturbations at recom- (Vol. 2), The University of Chicago Press, Chicago 
bination epochs. Therefore, though a better understand- 
ing of the physics underlying the resonant mechanism is 
needed, this point shall be thoroughly treated in a forth- 
coming paper of ours. 

Third: One must note that it will be interesting to 
study what happens for the case of a 3-component system 
like the dark matter mixed models. We suspect that in 
such a case more than one structure formation resonance 
occurs. This is a question which we intend to address in 
the future. 

Finally, one should also mention that although reso- 
nances between periodic coupled oscillators are common 
in all branches of physics, this type of resonance occurring 
in (non periodic) collapsing modes are, to our knowledge, 
yet unheard of, and we suspect that in other branches of 
physics such as in plasmas (where one has a 2-component 
fluid whose components are coupled by the electromag- 
netic interaction) similar resonances may perhaps occur, 
giving rise to instabilities. We think that such an open 
question is also worth investigating. 
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Fig. 1. Dependence on scale (in Solar Masses), of the eigen- 
value Al and imaginary part of A3. Note the two relevant mass 
scales for Ai and only one for A3 which coincides with the sec- 
ond one for Ai. The vertical arrow indicates the total mass M^, 
inside a spherical region of diameter 



Fig. 2. Typical projected phase-space trajectories, represent- 
ing the general solution of the dynamical system, on the plane 
(0:4,0:2). Also plotted are the projections of center and sad- 
dle subspaces. The left hand orbit corresponds to the situation 
Aj3^ < /?3 A^^ and the right hand one to the inverse initial 
condition. (This particular plot corresponds to a mass scale of 
the perturbation of ~ lO^^M©) 



Fig. 7. Coefiicients Ci and C2 and their dependence on scale. 
Note that the behaviour of these mode amplitudes, although 

not similiar to the one of ("3 and (4, is also peculiar for a scale 
around Ic- The vertical arrow indicates the total mass Mc inside 
a spherical region of diameter Ic 



Fig. 3. Projection on the (2:4, 0:2) plane of the three subspaces 
described in text for I > l^. The long-dashed lines represent 
the projection of the saddle subspace, in the limits I = Z„ and 
I — > 00. The short-dashed ones are the projections of the center 
subspace for the same limits 



Fig. 4. Scale dependence of the angle between the projection of 
the saddle and center subspaces with the X4-axis, respectively 
9i and ^3. The vertical arrow indicates the total mass M^, 
inside a spherical region of diameter 1^^ . Note the very peculiar 
behaviour around 



Fig. 8. Dependence on scale (in Solar Masses), of the growing 

mode amplitude ("3, for distinct initial conditions described by 
difl^ercnt values of the parameter q. One can notice that for 
q > 1 the values of (^3 to the right of the resonant scale are 
higher than their values to the left of it. The opposite occurs 
for g < 1 



Fig. 5. Coefficient Ca and its dependence on scale. The vertical 

arrow indicates the total mass Mc inside a spherical region of 
diameter Ic. Note the clear change in behaviour of this coeffi- 
cient around this scale 



Fig. 6. Coefficient ("4 and its dependence on scale. See Eq. (39) 
and the discussion following it. The vertical arrow indicates the 
total mass Mc inside a spherical region of diameter Ic 



Fig. 9. Dependence of the scale M, (in Solar Masses), at which 
the resonant peak in the amplitude Cs occurs, plotted as a 
function of the mass of the particles which constitute hot 
dark matter 
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